Tachyacoustic Cosmology: An Alternative to Inflation 
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We consider an alternative to inflation for the generation of superhorizon perturbations in the 
universe in which the speed of sound is faster than the speed of light. We label such cosmologies, first 
proposed by Armendariz-Picon, tachyacoustic, and explicitly construct examples of non-canonical 
Lagrangians which have superluminal sound speed, but which are causally self-consistent. Such 
models possess two horizons, a Hubble horizon and an acoustic horizon, which have independent 
dynamics. Even in a decelerating (non-inflationary) background, a nearly scale-invariant spectrum 
of perturbations can be generated by quantum perturbations redshifted outside of a shrinking acous- 
tic horizon. The acoustic horizon can be large or even inflnite at early times, solving the cosmological 
horizon problem without inflation. These models do not, however, dynamically solve the cosmolog- 
ical flatness problem, which must be imposed as a boundary condition. Gravitational wave modes, 
which are produced by quantum fluctuations exiting the Hubble horizon, are not produced. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION: WHY INFLATION WORKS 

Inflationary cosmology [1] |3] is the most successful 
and widely accepted paradigm for understanding the very 
early universe. By all measures inflation is a compelling 
and scientifically useful theory, and makes quantitative 
predictions which have been strongly supported by recent 
data [4j i5j j6j . The two main hallmarks of inflationary cos- 
mology are solutions to the flatness and horizon problems 
of the standard Big Bang cosmology: why is the universe 
so close to geometrically flat, and how did the apparent 
acausal structure of the universe arise? "Acausal" more 
specifically means that the universe is approximately ho- 
mogeneous on scales larger than a Hubble length H^^, 
and in addition exhibits a spectrum of density perturba- 
tions which is correlated on scales larger than a Hubble 
length. Such superhorizon correlations are generated in 
inflation by accelerated expansion, a/a > 0, where a(t) 
is the cosmological scale factor, which means that the 
comoving Hubble length dn — {aH)~^ shrinks with the 
expansion of the universe, 



din ( 



{aH)~^ < 0. 



(1) 



Therefore, quantum perturbations, which have constant 
wavelength in comoving units, are smaller than the Hub- 
ble length at early times, and are redshifted to larger 
than the Hubble length at late times, where they are 
"frozen" as classical perturbations. Furthermore, as long 
as the Hubble constant H is slowly varying with time, 
the perturbations generated in inflation are nearly scale- 
invariant, consistent with observation. Furthermore, the 



*dbessadaSbuf f alo . edu, 
t whkinneySbuf f alo . edu, 
* dsTTObuf f alo . edu, 
s jwangSniagara. edu 



solution to the horizon problem and the flatness problem 
are linked in inflation via a conservation law. 



\n-i\ 
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= 0. 



(2) 



Through this conservation law, a universe with shrinking 
comoving horizon size is identical to a universe which is 
evolving toward flatness. 
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- 11 < 0. 



(3) 



Inflation therefore solves the horizon and flatness prob- 
lems of the standard Big Bang with a single mechanism: 
accelerated expansion. 

However, inflation is not the only way to accomplish 
this goal, as can be seen from the fact that the acceler- 
ation a appears nowhere in the conservation law (|2|. To 
solve both the horizon and flatness problems, it is suffi- 
cient to have a shrinking comoving Hubble length. One 
way to do this is accelerated expansion, but another is 
to have a collapsing universe, H ~ (d/a) < 0. A collaps- 
ing, matter- or radiation-dominated universe also has a 
shrinking comoving Hubble length, which will generate 
perturbations in a manner similar to inflation. This is 
the mechanism used by the Ekpyrotic scenario ^ to con- 
struct a cosmology consistent with observations. It is also 
possible to decouple the horizon and flatness problems, 
for example in theories with a varying speed of light, so 
that the causal horizon is much larger than the Hubble 
length [8 . Such theories can in principle solve the horizon 
problem, but not the flatness problem, since the conser- 
vation law (2]) is violated. It is also possible to solve the 
horizon problem by a universe which is much older than 
a Hubble time as in string gas cosmology [TT] or island 
cosmology [5J [TD], or by the inclusion of extra dimen- 
sions [12, 13]. However, it has been argued that inflation 
and ekpyrosis are the only mechanisms for generating a 
scale-invariant spectrum of perturbations [71 114j. 
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In this paper, we discuss a method of solving the cos- 
mological horizon problem and seeding scale-invariant 
primordial perturbations in a cosmology with deceler- 
ating expansion and a corresponding growing comoving 
Hubble horizon. The key to implementing such a model 
is the fact that curvature perturbations are not gener- 
ated at the Hubble horizon, but at the acoustic hori- 
zon determined by the speed of sound of a scalar field. 
For canonical field theories, the two are identical, but 
for non-canonical field theories, they are not. If one 
has a decaying, superluminal sound speed, curvature per- 
turbations can be generated outside the Hubble horizon 
without inflation. We propose the term tachyacoustic 
for such cosmologies, which are closely related to varying 
speed of light theories. This idea has some history: such 
cosmologies were first proposed by Armendariz-Picon in 
the context of modified dispersion relations |15j , and the 
generation of perturbations in such cosmologies was fur- 
ther considered by Piao [16^ . The idea re-emerged in the 
context of varying speed of light theories by Magueijo 
[17j . and non-canonical Lagrangians by Magueijo [TS] 
and Piao [12]. In this paper, we outline a general ap- 
proach to such cosmologies based on the generalization of 
the infiationary fiow formalism [20j introduced by Bean, 
et al. for the case of arbitrary Lagrangians [21]. We 
find that there is a class of Lagrangians with the neces- 
sary properties for tachyacoustic cosmologies, and discuss 
two interesting examples. We find that it is straightfor- 
ward to generate nearly scale-invariant perturbations for 
these Lagrangians, and show that they have the prop- 
erty of reducing to instanton-like solutions with infinite 
sound speed on the initial-time boundary of the space- 
time. We speculate that this property may allow a 
self-consistent description of tachyacoustic cosmologies 
within a Wheeler-deWitt description of quantum cosmol- 
ogy. Finally, we show that such models are causally self- 
consistent, and argue that they form a viable class of 
alternatives to inflation. 



II. TACHYACOUSTIC COSMOLOGY 

In this paper, we consider a way of generating scale- 
invariant superhorizon cosmological perturbations based 
on non-canonical scalar field Lagrangians with a speed of 
sound faster than the speed of Hght, cg > 1. If the uni- 
verse is dominated by a scalar field with speed of sound 
cs , the relevant horizon for the generation of density per- 
turbations is not the Hubble horizon dn — {aH)^^ but 
the acoustic horizon. 



Dh 



aH' 



(4) 



Mode freezing at the acoustic horizon is well-known in 
non-canonical inflation models, for example k-Inflation 
[22] and DBI inflation [23]. In non-canonical inflation 
models, the Hubble horizon and the acoustic horizon are 
both shrinking in comoving units, resulting in the gen- 
eration of density perturbations at the acoustic horizon 



and gravitational wave perturbations at the Hubble hori- 
zon [24j . However, the comoving Hubble horizon need 
not be shrinking to generate curvature perturbations: all 
that is required is that the acoustic horizon be shrink- 
ing, dDn /dlna < 0. In this case, if curvature perturba- 
tions are to be generated on scales larger than the Hub- 
ble horizon, it is necessary that the acoustic horizon be 
larger than the Hubble horizon, which requires a speed 
of sound greater than the speed of light. Such theories 
were studied recently by Babichev et al. [25j |26], who 
showed that k-essence theories with > 1 are causally 
self-consistent (see Appendix [A| , and can be mapped to 
bimetric theories with two "light cones" , one given by 
the Hubble horizon, and the other given by the acous- 
tic horizon, which can be larger than the Hubble hori- 
zon without the presence of closed timelike loops. This 
opens the possibility that one can construct a deceler- 
ating cosmology which nonetheless generates perturba- 
tions on super-Hubble scales via a superluminal, shrink- 
ing acoustic cone. 

To explicitly construct such a model, consider a Dirac- 
Born-Infeld (DBI) Lagrangian, 



c {cf>, X) = -r' wyi-f (0) cb^ + r' w -v{<f>), 

(5) 

where we take the metric to have negative signature 
Vfif = diag.(l, — 1, — 1, — 1), so that the kinetic term for 
the field is positive. 



02 EE g^''d^<f>d,^ > 0. 



(6) 



Such Lagrangians often arise in string theory, for example 
in the case of DBI inflation. For the moment, we will not 
attempt to make a connection with string theory, but will 
take the Lagrangian (|5]) as a phenomenological ansatz. 
Take the background spacetime to be a flat, Friedmann- 
Robertson- Walker (FRW) space, = a'^'q^i,, so that 
the scale factor evolves as 



a oc exp 



Hdt 



-N 



(7) 



where we deflne the number of e-folds A^ as^ 

N=- I Hdt. (8) 



There is a class of exact solutions [271 128] to the equation 
of motion for the field characterized by two dimension- 
less fiow parameters e and s, where 



and 



1 dH 



1 dcs , 

s — — = const. 

Cs dN 



(9) 



(10) 



^ We use the usual convention that N —* oo corresponds to early 
time, and N —> — oo corresponds to late time. 
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The parameter e has its usual interpretation in terms of 
the equation of state of the scalar field, 



P 



(11) 



For e = const., the scale factor evolves as a power-law, 
a cx t^/'^, so that the expansion is accelerating {i.e. infia- 
tion) for e < 1. The speed of sound evolves as 



and the Hubble parameter evolves as 



H = - (X e 
a 



(12) 



(13) 



The parameter e is a positive-definite quantity for p > 
—p, so that the Hubble constant always decreases with 
expansion. In contrast, the parameter s can take either 
sign, with s > corresponding to a sound speed which 
increases with expansion, and s < corresponding to an 
decreasing sound speed. (See Ref. [55] for a detailed 
derivation of this solution.) The important dynamics for 
the generation of perturbations is the time evolution of 
the corresponding horizons in comoving units. The co- 
moving Hubble horizon evolves as 



dn oc (aH) 



oc e(i-^)^ (X T, 



(14) 



where r is the conformal time. The Hubble horizon is 
shrinking in comoving units for e < 1, which is identical 
to accelerated expansion, and is the usual condition for 
inflation. The acoustic horizon behaves as 



oc — cx e 
an 



(1- 



(15) 



Therefore the condition for a shrinking acoustic horizon, 
1 — e — s > 0, is not identical to accelerated expansion. 
For e > 1 and s < 1 — e, the expansion is non- inflationary, 
the Hubble horizon is growing in comoving units, and the 
acoustic horizon is shrinking. The initial singularity is at 
T = 0, and we see immediately that for the tachyacoustic 
solution, the speed of sound in the scalar fleld is infi- 
nite at the initial singularity, and the acoustic horizon 
is likewise inflnite in size. Therefore, such a cosmology 
presents no "horizon problem" in the usual sense, since 
even a spatially inflnite spacetime is causally connected 
on the initial-time boundary. Furthermore, unlike in the 
case of inflation, there is no period of reheating neces- 
sary, since the cosmological evolution can be radiation- 
dominated throughout and the cosmic temperature is not 
driven exponentially to zero. 

In the next section, we use the the generalized flow 
function approach of Bean, et al. [2T] to construct a 
class of Lagrangians with solutions of the type outlined 
above, with constant flow parameters. In these solutions, 
the scale-factor evolves as a power-law in time and the 
equation of motion for curvature perturbations can be 



solved exactly, which we discuss in Sec. VI 



III. FLOW HIERARCHY FOR GENERAL 
K-ESSENCE MODELS 

We now generalize the discussion in the last section 
to an arbitrary k-essence model. Consider a general 
Lagrangian of the form £ = £ [X, (f>] , where 2X = 
g^^df^(t)d^(l) is the canonical kinetic term {X > accord- 
ing to our choice of the metric signature). The energy 
density p and pressure p are given by 



P 
P 



c{x,4>), 

2XCx - C. 



The speed of sound is given by 



4 



Px 

PX 



1 + 2X 



XX 



X 



(16) 
(17) 



(18) 



where the subscript "X" indicates a derivative with re- 
spect to the kinetic term. Throughout this section, unless 
otherwise stated, we will follow closely Bean et al. |21j . 
We deflne the flrst three flow parameters as derivatives 
with respect to the number of e- folds, dN = —Hdt:^. 



1 dH 
HdN' 

1 dcs 

cs dN' 

1 dCx 
Cx dN ■ 



(19) 
(20) 
(21) 



The Friedmann equation can be written in terms of the 
reduced Planck mass Mp = l/%/8^ 



1 



1 



p = 



Mil 



{2XCx - C) , 



and the continuity equation is 



p = 2HH = -3H ip + p) = -GHXCx ■ 



(22) 



(23) 



For monotonic fleld evolution, the fleld value (f> can be 
used as a "clock" , and all other quantities expressed as 
functions of (j), for example X = X {(j>), C = £[X {(j>) , (j)], 
and so on. We consider the homogeneous case, so that 
(/) = V2X. Using 



d 
dt 



(24) 



we can re-write the Friedmann and continuity equations 
as the Hamilton Jacobi equations. 



X 



-c. 



X 



(25) 
(26) 



^ The parameters s and s correspond to the parameters ft and k 
in Bean, et al. |21| 
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where a prime denotes a derivative with respect to the 
field 4>. The number of e-folds dN can similarly be re- 
written in terms of d(j) by: 



dN 



'Hdt 



H 

/2X 



2M|, 



del) 

H' (0) 



d(j). 



(27) 
(28) 



IV. COSMOLOGICAL SOLUTIONS FOR 
CONSTANT FLOW PARAMETERS 

The simplest way to solve the flow equations derived in 
the preceding section is to take all of the flow parameters 
to be constant, 



(k_ 
dN 



ds 
dN 



ds 
dN 



d^ 
dN 



da 
dN 



d^ 
dN 



0. 



(34) 



The flow parameters e, s, and s ([19]) can therefore be Then, from ( [T9pT| we easily find the following relations: 



written as derivatives with respect to the field (j) as: 



2Ml 

2Mj, H' (0) c's (0) 
Cx H{4,) C5(0)' 
2APp H' (0) C'x 
Cx H{<P) Cx ■ 



(29) 
(30) 
(31) 



Taking successive derivatives d/dN with respect to the 
number of e-folds yields an infinite hierarchy of flow equa- 
tions [Mun], 



de^ 

dN 
dfj 

dN 
ds 

dN 

ds 

dN 
d^ 

dN 
d^ 

dN 

d^(3 

dN 



(2e - 27] + S) , 
-77 (e + S) + 

{e — 7] + s + s) + ep, 

-S{e- ii + 2S) + e^P, 

~''a[{e-l){e~fi)+lS + s]+''+^a, 
J(3[{£-l){e-fi) + {£+l)S]+'+'(3, (32) 



where the higher-order flow parameters are defined as 
follows, where ^ = 1, . . . cxd is an integer parameter: 



77(0) 

'A(0) 
^a(0) 



^ 2MlH"{4>) 
Cx H{<t>) ■ 
f2Ml\^ fH'{(j)) 



f2Ml 
/2Af|, 



X 



X 



H{<j>) 



1 



d'+- 



H {(t>) d<j>'^ 
1 d^+i 



1 d^+i 

-Cx. (33) 



C 



X 



Solutions to this infinite hierarchy of flow equations are 
equivalent to solutions of the scalar field equation of mo- 
tion. In the next section, we specialize to the case where 
the flow parameters are constant, which results in an ex- 
actly solvable system. 



H (X e 
cs oc e 



eN 
-sN 



Cx oc e''\ (35) 
The first two are identical to the DBI case, Eqs. 



(121 



and (13 1, but in the fully general case Cx evolves inde- 



pendently of Cs- It is straightforward to verify that the 
full flow hierarchy ( 33 1 reduces to an exactly solvable set 



of algebraic equations, with the higher-order parameters 
expressed in terms of e, s, and s. We can use the rela- 
tions ( 29|30|31 1 to solve fo r H {(j)), cs (0), and Cx {4>) as 
follows: from Eqs. ( |29|31[ ), we have 



2Ml 



H' 
1[ 



2e^^ = const. 
Cx 



Cx \n J L,x 
We then have a differential equation for Cx 
C'x 



Mp^/2e 



const . , 



with solution 



Ix (<^) = ^ 



e fMp 



(36) 



(37) 



(38) 



where the integration constant has been absorbed into a 
field redefinition. From Eq. (35 1, the field then evolves 
as 



sN 



(39) 



so that the direction of the field evolution depends on the 
sign of s. 



'^4^-~-dN. 
6 2 



Equation ( 29 1 then reduces to 
ll 



4e^ 



with solution 



H cx 



,±2£/s 



(40) 



(41) 



(42) 



The sign ambiguity can be resolved by requiring that the 
universe be expanding, dH/dN > 0, so that 



H cx 



cx e'^. 



(43) 
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Finally, we solve for the speed of sound using Eq. ( 30 1 
which reduces to 



2s 

— — = const. 



with solution 



(44) 



(45) 



Since our choice of A'^ = corresponds to an arbitrary 
renormalization of the scale factor a oc , we can with- 
out loss of generality define C5 = 1 at = 0, so that the 
general solution for the background evolution is given by 



cs (0) = 

where the field evolves as 



Mf 



-2e/s 



2s/s 



-sN/2 



(46) 



(47) 



(48) 



(49) 



We can derive the time dependence of the scale factor 
using the Hamilton- Jacobi equation (251, 



= -H{^)^^V2X, 
so that the kinetic term can be written as 

Integrating expression ( [50| gives 

et 

so that the scale factor evolves as a power-law in time, 
consistent with the relation ^ between e and the equa- 
tion of state w — p/p, 



(50) 



(51) 



(52) 



.(t)ocii/^-t2/3(i+»)^ 



(53) 



Radiation-dominated evolution therefore corresponds to 
e = 2, and matter-dominated evolution corresponds to 
e = 3/2. Infiation corresponds to e < 1. The comoving 
Hubble horizon evolves proportional the to the conformal 
time, 



1 (X e(i-^)^ cx T, 



dn oc {aH) 
and the acoustic horizon evolves as 



aH 







(54) 



(55) 



identically to the DBI case discussed in Sec. |IT] For e > 
1, the acoustic horizon is shrinking in comoving units 
for s < 1 — e. Note that this behavior is independent 
of the parameter s, which determines the form of the 
Lagrangian, as we discuss in the next section. 



V. RECONSTRUCTING THE ACTION 

In the past two sections we have solved the fiow hier- 
archy for a model characterized by constant flow param- 
eters, which allowed us to solve for H{(j)), cs ((/)), and 
Cx {4')j only the derivative of the Lagrangian with re- 
spect to the kinetic term X is determined. Therefore 
this solution corresponds not to a single action but a 
class of actions. In this section we derive a general equa- 
tion for Lagrangians in this class, and discuss two specific 
examples. 

From Eqs. (38) and (48 1, we see that the speed of 



sound Cs can be written in terms of Cx 



4 



q—\q — 2sIs 



1 + 2X 



C 



XX 



C 



X 



where we have used Eq. (18 1, and defined 



C . 



(56) 



(57) 



The result is a differential equation for the function 
C{X,<j,): 



2XCxx + Cx- CC\ = 0, 



where we have defined 



1 



2s 

s 



(58) 



(59) 



Therefore, by specifying a relationship between the pa- 
rameters s and s, we can construct a Lagrangian as the 
solution to the differential equation (58). For example, a 
canonical Lagrangian with speed of sound C5 = const. = 
1 is just the case s = 0, so that n = 1 and C — \, and 
Eq. (58) becomes 



-XX 



0, 



with general solution 



£ = /(0)x-y(0) 



(60) 



(61) 



Here / (i/i) and V (0) are free functions which arise from 
integration of the second-order equation ( 58 1 . The func- 



tion / (0) can be eliminated by a field redefinition d^p = 
f {(j))d(l), resulting in a manifestly canonical Lagrangian 
for if, as we would expect from setting cs = 1. A canon- 
ical Lagrangian can support infiationary solutions, but 
not tachyacoustic solutions, and is therefore not of inter- 
est here. However, other choices of n do yield tachya- 
coustic solutions, and we focus on two such choices: 

1. 71 = 0: A Cuscuton-like model. 

2. 71 = 3: A DBI model. 

We discuss each case separately below. 
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n = 0: A Cuscuton-like model 



The Hamilton- Jacobi Equation (25) becomes: 



The case n = corresponds to s = —2s in (59 1, with 
solution 



(62) 



£(X,0) = 2/(0) ^X + CX^V{ci>) 



This Lagrangian is similar to a "cuscuton" Lagrangian 
[36j . with the addition of a term proportional to X. 
Unlike the original Cuscuton model, which represents a 
causal field with infinite speed of sound, the solution ob- 
tained here is valid for the general case, in which the 
speed of sound can be finite. A similar cuscuton-like La- 
grangian was considered in Rcf. |19j . 

As in the canonical case, the functions / {(f) and V (0) 
are free functions resulting from integrating Eq. ( 58 1 . 



Unlike the canonical case, however, neither can be re- 
moved by a field redefinition. However, both functions 
are fully determined by our choice of solution with e, s. 



and s constant. Differentiating Eq. (62) with respect to 
X gives 



2e f Mp 



(63) 



where the right hand side is the solution ( 46 1 . Then 

2 " 



2e 



Mp 
Mp 



1 



(64) 



where for 2s — — s, the expression (48) for the speed of 
sound becomes 



cs (0) 



(65) 



The Lagrangian ( 62 1 can then be written as 

2 



[24 ((/.)-!] -1^(0). (66) 



The Hubble parameter ( 47 1 is given by 

e/s 



ff(0) = iJo 



"0 



(67) 



and we can then express the kinetic term as a function 



of <j) using Eq. (51 ) 



^ ^ 2 ^ 2 c2(0) 



(68) 



The Lagrangian ( |62| can then be written entirely as a 
function of the field (h. 



■^MIh'^ = 2Mj,eH^ - C 

= y(0) + ^2 — ^ 



and we have an expression for the potential V {(j)), 



(70) 



(71) 



The Hubble parameter H (0) and the speed of sound 
Cs {(j>) are given by Eqs. (|67[) and (65), respectively. For 



/)o ^ 1, the speed of sound is much greater than the 
speed of light, C5 ^ 1, and the potential is approximately 



V (0) ~ 3M|,7J2 (^) = 3M|,i?2 



2e/s 



(72) 



which can be recognized as a slow-roll-like solution dom- 
inated by the potential ~ V^/3M^. For s < 0, the 
field is rolling away from the origin, and for s < 1 — e the 
comoving acoustic horizon is shrinking and the solution 
is t achy acoustic. 



B. n = 3: The DBI model 



The case n = 3, corresponds to s = s; then, from (30) 
and (31), we find £x = Cg^- Eq. (59) is then 



(73) 



so that we can take C ~ 1 without loss of generality. 
Therefore, the Lagrangian assumes the well-known DBI 
form, 



£ {X, 4>) = -r' (0) V^-f{4>)x + /-I ((/.) -!/(</.). 

(74) 

The DBI model with constant flow parameters is exten- 
sively discussed in Ref. [28 , and the reader is referred 
this paper for further details. For e and s constant, the 
functions V and / are fully determined and are given by 



y(0) = 
/(0) = 



3Mj,H^ (0) 
/ 1 



2e 

y 

1-4(0) 



1 



1 + cs (0) 



\2Ml,eJ H^^)cs{^)' 



(75) 



The Hubble parameter and speed of sound are given by: 

-2e/s 



(76) 



and 



C = Mf,eH^ (0) 



1 



4 (0) 



V (0) . (69) 



cs (0) 



(77) 
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DBI Lagrangians allow for either inflationary or tachya- 
coustic evolution [18j . depending on the values of e and 
s. Note that for cg > 1, the function / is negative, which 
has consequences for embedding such a model in string 



theory, which we discuss in Sec. VII 

In this section, we have explicitly constructed La- 
grangians, including fully determined potentials, for 
which the flow parameters are constant and the back- 
ground evolution can be solved exactly. For suitable 
choices of the flow parameters, the evolution is tachya- 
coustic, i.e. with a growing comoving Hubble horizon and 
a shrinking comoving acoustic horizon. In the next sec- 
tion, we discuss the generation of curvature perturbations 
at the acoustic horizon and show that such perturbations 
are nearly scale-invariant, consistent with observation. 



where a prime denotes a derivative with respect to con- 
formal time, ds^ = a? (r) (dr^ — dx^) . It is easy to show 
that the variable z, defined by ( [sT] ) can be cast into the 
following form. 



aMpy/2e 
cs 



then, using 



d 
dr 



-aH- 



dN' 



(83) 



(84) 



we can evaluate the ratio z" / z in (82 1 in terms of the 
flow parameters ( 29][33| ; the result is 



(85) 



VI. COSMOLOGICAL PERTURBATIONS FOR 
CONSTANT FLOW PARAMETERS 

We can deal with cosmological perturbations in this 
general k-essence model with constant flow parameters 
in the same way as performed in |28j . Following the ap- 
proach of Garriga and Mukhanov [21] we start with the 
perturbed Einstein equations. 



d 
dt 



= $ 



a2(p + p) 

a{p + p) 
2Ml 



where <i> is the Bardeen potential and 
bation of the fleld 



(78) 



is the pertur- 



Equations ( 78 1 can be cast into a 



more convenient form by changing the perturbations <i> 
and Scf) to the new variables C, and ^ deflned by 



H 



(79) 



so that the perturbed Einstein equations ( 78 ) become 



• ^ a{p + p) 

C = — ^ v^^. 



(80) 



As usual, we introduce a new variable z and the gauge- 
invariant Mukhanov potential u as 



a(/9-|-p)i/2 



u = zC; 



(81) 



then, from (80) we derive the mode equation for u{t) oc 
Ufc(r) exp(ik • x), given by 



(csfc)' 



Uk = 0, 



(82) 



where 



F = 2 + 2e-3f)-3s+ ^S + 2e2 + 

2 4 



2ss 



+ fj'^ + 2e(s — s) + Sijs — —rjs — Af]e + 



(86) 



Next, it is convenient to change the conformal time, r, 
to the ratio of wavenumber to the sound horizon. 



csk 

then, conformal time derivatives switch to 



(87) 



and 

d^ 
dr^ 



d 

d^ 



(1 



-aH{l 



\2 2 

s) y 



dy'^ 



dy 



where 
G 



d 

dy 
(89) 



-s + e(2s + s) + s(2s + s) + 2e^ - 2ef] - sfj 



- e (^a) . (90) 

It is important to stress that the functions F and G 
derived above hold in general; they reduce to the well 
known expressions in the DBI limit |28j . which, in this 
case, s — s and ^a = -^f] = p. Substituting (85l and (89l 
into the mode equation ( 82 1 , we flnd 



dy'^ 



Gy 



duk 
dy 



[y 



F]uk^O, (91) 



which is an exact equation, without any assumption of 
slow-roll. 

In the case where the flow parameters are constant, 
we can use the differential equations ( 32 1 to reduce the 



number of independent parameters. We have 
1 



1] 



(2e + S), 



8 



'a 



2e 
2e ' 



(2s + ~s) , 



(92) 



then, substituting these values into expressions (86 1 and 
(90 1, we find, respectively. 



9 n 3 

F = 2 - e - 3s + -s^ - -ss + es 



G = s(-l + e + s). 



(93) 



(94) 



It is important to notice that F is different from the cor- 
responding expression found in the DBI case [28] , since 
the gauge-dependent s comes into play. However, G is 
identical to its DBI analog, and it is expected since basi- 
cally it comes from the change of variables r — > y, which 
depends solely on the parameters cs and H, and not on 
Cx- For constant flow parameters we can solve equation 



( 91 1 exactly, and the solutions are given by 



V 



(1) 



1-e-s 



1-e-s 



(95) 



where Ci and C2 are constants, and H^p , H^'^^ are Hankel 
functions of first and second kind, respectively. The order 
v of the Hankel function is given by 



9 - 6e - 12s + 9s^ - 3sS + 4es - 2s 
4(l-e-s)2 



2 , ,2 



next, using (35), (|84[) and (88 1 we find that 



Cs oc y 



s/(e+s-l). 



(96) 



(97) 



then, imposing the Bunch-Davies vacuum C2 = 0, and 
normalizing the mode amplitudes by means of the canon- 
ical quantization condition 



»duk 
^ dy 



dul 



Uk 



dy Csk{l 



=s)' 




y 



y 



(98) 



(99) 



1 — e— s "V^l — e — s^ 
which differs from the DBI case only in the order of the 



Hankel function ( 96 1 . In the small wavelength limit y — > 
00 the early-time behavior of Uk will be identical to DBI 
P5] for constant flow parameters 



1 



piy/il-e-s) 



(100) 



whereas in the late-time behavior y ^ the mode func- 
tion behaves as 



,,1/2-1^ 



From (101) we can derive the expression for the scalar 



spectral index . Using the definition of the power spec- 
trum of curvature perturbations 



27r2 



(102) 



and substituting expressions (83) and (101) into (102) 
we find 



n 



(103) 



87r2M|, cse 

at horizon crossing, where Jiy) is a constant given by 

2-3/2^(1 -.-s)-V2. (104) 

then, from the definition of the scalar spectral index 



1 



din A: 



(105) 



and using 



1 



d 

SI dN' 



dink \1 — e 
we see that the spectral index assumes the form 

Tic = 1 



(106) 



1-e-s' 



(107) 



which does not depend on the gauge-dependent parame- 
ter s, and is identical to its DBI analog. This is expected 
since the power spectrum evaluated at the horizon cross- 
ing, equation (103), depends solely on H and cs, whose 



derivatives with respect to N are related to the gauge- 
invariant flow parameters e and s. The scale-invariant 
limit is s = — 2e. 



VII. CONCLUSIONS 

In this paper we have demonstrated that accelerated 
expansion or a collapsing universe are not the only ways 
to dynamically generate a scale-invariant spectrum of 
superhorizon curvature perturbations. There is a third 
way: a superluminal acoustic cone which is shrinking in 
comoving coordinates. Curvature perturbations gener- 
ated at the acoustic horizon are familiar from inflation- 
ary scenarios based on non-canonical Lagrangians such 
as /c-inflation and DBI inflation. Such non-canonical La- 
grangians arise naturally in string theory. However, in 
these scenarios, both the Hubble horizon and the acoustic 
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horizon are shrinking in comoving units, and the acoustic 
horizon is typically smaller than the Hubble horizon, i.e. 
C5 < 1. It is natural to ask whether tachyacoustic mod- 
els have a similar, natural stringy embedding, especially 
since the DBI action Q naturally admits tachyacous- 
tic solutions. Such an embedding is nontrivial, however, 
since the frequently considered case of a 3-1-1 dimension 
d-brane evolving in a higher-dimensional throat is ill- 
defined in the cs > 1 limit. To see this, consider the 
full ten-dimensional metric of throat plus brane [29j , 

dsfo = (r) dsj + h-^ (r) {dr^ + r^ds\^) . (108) 

The field cf) is simply related to the coordinate in the 
throat r as (j) — \/T^r, where the brane tension de- 
pends on the string scale and the string coupling gs 
as [35] 



(109) 



The Lagrangian for the field can be shown to be of the 
DBI form ([s]), where the inverse brane tension /(</)) is 
given in terms of the warp factor /i ((/>) by 



where we have s < — 3 for a shrinking comoving acoustic 
cone. For s < —4, the spectral index is blue, n > 1, 
which is ruled out by observation. The WMAP 2a limit 
n = 0.96 ±0.026 [5] corresponds to s = [-3.814, -3.959]. 
Since the Hubble horizon is growing in comoving units, 
no gravitational wave modes are produced. 

Tachyacoustic models are not a fully convincing alter- 
native to inflation, since they solve only the horizon prob- 
lem and not the flatness problem, and inflation solves 
both at once. However, inflation has initial conditions 
problems of its own, in particular the fact that the initial 
inflationary "patch" must be larger than a horizon size 
for inflation to start |37j . Furthermore it has been shown 
that inflationary spacetimes are in general geodesically 
past- incomplete [38'. The initial conditions for tachya- 
coustic cosmology are quite different than those for infla- 
tion due to the presence of a true "Big Bang" singularity 
at zero time. However, in this limit, the sound speed is 
infinite and the tachyacoustic solution approaches an in- 
stanton. To see this, examine the form of the DBI field 
Lagrangian (|5| near the t — boundary of a tachya- 
coustic spacetime. From Eq. (12), the C5 — > 00 limit 
corresponds to — > 00 and / (cf)) 0^ — > —00, so that 



/(0) 



(110) 



The problem is immediately evident: super luminal prop- 
agation C5 > 1 requires / < 0, so that the factor (0) 
appearing in the metric (108) is imaginary, and the met- 



ric is ill-defined. Therefore, although the DBI action it- 
self admits tachyacoustic solutions, this limit does not 
correspond to a well-defined string solution. It is not 
clear whether or not string manifolds exist which self- 
consistently admit solutions with C5 > 1. 

We calculate the scalar spectral index of perturbations 
for tachyacoustic solutions, and find 



1 



2e 



1 



(111) 



Unlike inflationary models, radiation-dominated tachya- 
coustic models do not require a period of explosive en- 
tropy production to transition to a "hot" Big Bang cos- 
mology. The early universe must be scalar-fleld domi- 
nated, but the temperature of the universe is not driven 
exponentially to zero, since the scalar has a radiation 
equation of state at all times, and entropy density is con- 
served (for any radiation component with density p-y, the 
ratio Ptji/py = const). The scalar field must eventually 
decay to Standard Model degrees of freedom, but as long 
as this happens before primordial nucleosynthesis, the 
model will match observations. A slow or late decay of 
(j> into other degrees of freedom would also suppress the 
production of unwanted relics such as monopoles or grav- 
itinos. For radiation-dominated tachyacoustic expansion 
with e = 2, the spectral index is 



n= 1 



(112) 



C 



y(0) 



where 



(113) 



(114) 



From Eq. ( 75 1 , the asymptotic behavior of V (0) and 
/ {(j}) are 



V{<f>) 



3Mj.H^ cx 
1 cs 
2Mle H 



-it/s 



cx 



2(l+2e/s) 



(115) 



The scale-invariant limit s = — 2e is especially interesting, 
since 



(116) 



^ p = const., 

and the Lagrangian takes the form 

C^p^-V (0) , 



(117) 



where V {(p) oc (fP' . This can be identified as exactly the 
"cuscuton" Lagrangian, suggested by Afshordi, et al. as a 
candidate for Dark Energy |36ll40JST]. Similarly, the n = 
solution considered in Sec. |V] approaches a cuscuton 
on the initial boundary surface. The cuscuton is a non- 
dynamical, instanton-like solution with infinite speed of 



sound. Consider the action for the Lagrangian (117) 



Sri 



p^-V{cj,) 



p 



dt(j) 



d^x^ 



5- / d'^xV^VicP) 
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dcjjY. ((/)) - 



-gv (0), 



(118) 



where S ((f)) is the volume of a constant-0 hypersurface 
in the spacetime. The classical solutions to the cuscuton 
action are constant mean curvature hypersurfaces, anal- 
ogous to soap bubbles ^36j . It is interesting to speculate 
that this property of the cuscuton action may provide a 
self-consistent cosmological boundary condition, or (even 
more speculatively) be useful as a solution to the cosmo- 
logical flatness problem. A full analysis, however, would 
require inclusion of the gravitational action and solution 
in a Wheeler-Dc Witt framework, or perhaps an embed- 
ding of the model in string theory or an alternate gravity 
theory such as Horava-Lifshitz [42l |43] . This is the sub- 
ject of future work. 



APPENDIX A: STABLE CAUSALITY 

Since tachyacoustic cosmology deals with superluminal 
propagation of perturbations, it is important to address 
the issue of causality in this model. Babichev et. al. [26 
have discussed the conditions that must be fulfilled by a 
general k-essence model with superluminal propagation 
in order to avoid causal paradoxes (i.e., the presence of 
closed causal curves - CCC). In this appendix we outline 
the main ideas of this work and apply to our tachyacous- 
tic model. 

To begin with let us introduce some key definitions 
[44j . Let g^u be a metric with Lorentzian signature de- 
fined on a given manifold M.. Given a point p € Al, 
let be a timelike vector at p\ then, from this timelike 
vector we construct a second metric, g^v, related to the 
background metric g^i, by 



(Al) 



which clearly has a Lorentzian signature. Then, the 
spacetime {M,gp,^) is said to be stably causal if there 
is a continuous timelike vector field t^ such that the 
spacetime (7U,5^i/) possesses no closed timelike curves. 
The following theorem (8.2.2. in [H]) establishes the 
necessary and sufficient conditions for a spacetime to be 
stably causal: 

A spacetime {Jv[,g^^) stably causal if and only if 
there exists a differentiable function f on Ai such that 
"S/^f is a past directed timelike vector field. 

We can apply this theorem to k-essence models as 
follows [26]. First, we must find the analog of the 
induced metric (All for the case of k-essence models, 
which can be obtained by means of the equation of 
motion for a scalar field described by a Lagrangian 



£(X,0), 



2XCx^ - = 0, (A2) 
where G^"^, called "effective" or "acoustic" metric is given 

by 



It is convenient to use the metric |2S1 



G^" 



'-x 



(A3) 



(A4) 



fines the same causal structure. The inverse metric G^J 



which is conformally equivalent to G'^'^, and hence, de- 
fines the sai 
is given by 



^-1 ^ 

cs 



gfiv — Cg ■ 



-V^</)Vi.<; 



(A5) 



notice that it has the same form of (All, since is 



timelike. Using this definition, we can now apply the 
theorem stated above and check the stable causality of 
k-essence models. Let t be time coordinate with respect 
to the background metric (which is everywhere future 
directed), which we take to b e FRW. Since Vf^tV^t = 1, 
we have, using (A3 1 and (A4|, 



cs 



C 



X 



C 



XX ; 2 



C 



X 



(A6) 



then, since for a homogeneous scalar field holds cjp' — 2A, 
we have, from (18) and (A6l that 



csC 



>0, 



X 



(A7) 



provided the Null Energy Condition (NEC) is satisfied, 
that is, Cx > 0. Therefore, t plays a role of global time 
for both spacetimes {A4,gij,^) and (A4,G~^), and then 
the conditions of the theorem are fulfilled. Then, there 
is no CCC in superluminal k-essence models built from 
homogeneous scalar fields on a FRW background. Since 
this is exactly the case of the models introduced in this 
paper, we conclude that there are no causal paradoxes in 
tachyacoustic cosmology. 
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